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Abstract 

It has been argued recently that transverse asymmetries that are expected to be shielded from 
the presence of the S-wave Kir pairs originating from the decay of a scalar Kq meson, are 
indeed affected by this pollution due to the impossibility to extract cleanly the normalization 
for these observables. In this short note we show how using folded distributions, which is 
nowadays the preferred method to obtain the information from the 4-body decay mode 
4^ • B — > K*(—t Kn)l + l~ , one can easily bypass this problem and extract the clean observables 

-Pi, 2. 3 and also P4 5 6 in a way completely free from this pollution including all lepton mass 
corrections. We also show that in case one insists in using uniangular distributions to extract 
these observables it is possible to reduce this pollution to just lepton mass suppressed terms. 
On the contrary, the Si observables, that are by definition normalized by the full differential 
decay distribution, will indeed suffer from this pollution via their normalization. Finally, 
we also present a procedure to minimize the error associated to neglecting lepton mass 
corrections in the distribution defining a massless-improved limit. 
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The angular distribution B — > K*(—> Ktt)1 + 1~ is nowadays, but also in the near future, 
one of the central players in our search for New Physics in rare B decays. A wide variety of 
experimental analysis have been presented [U O El |4] focusing on different observables related 
to this rich mode. Simultaneously, there has been a huge effort from the theory community 
to provide a complete and accurate description of this 4-body decay mode [5]-[16]. Given the 
relevance of this mode in providing new constraints |17j-|28j it is of utmost importance to control 
and isolate any possible source of hadronic pollution that can spoil a clear signal of the theory 
that lies beyond the SM. 

In a recent interesting paper [16] one possible source of pollution of the angular distribution 
B — > K*(—> Kir)l + l~ due to events coming from the distribution B —> Kq{-^ Kit)1 + 1~ , where 
Kq is a scalar meson resonance, was analyzed. There it was presented the full distribution of the 
combined channels B — > K*(—t Kir)l + l~ and B — > Kq(—> Kit)1 + 1~ and the focus was on the 

(2) 

impact that this pollution could have on the extraction of the transverse asymmetries A T [6], 
|14j and j4^ m ^[I3]. In |16| it was argued that even if one should expect that the transverse 
asymmetries are unaffected by this S-wave contamination, they are afflicted from this disease 
via their normalization. They conclude that if the transverse asymmetries (defined with an 
unusual normalization) are measured in the particular way they propose, they are afflicted by 
this pollution at the level of less than 10% in all c/ 2 -range except for q 2 = 2 GeV 2 where the 
pollution can be as large as 23%. 

In this short note we present a simple procedure to bypass this pollution and extract the 
transverse asymmetries defined as in [15] in a way completely free from any S-wave contamination 
including also all lepton mass corrections. First, we show using folded distributions, that an 
exact (m,£ ^ 0) and clean (no S-wave pollution) extraction of all the P, observables [15] is 
perfectly possible. Second, we argue that even using uniangular distributions a strategy can be 
designed to reduce the impact of the companion decay on the transverse asymmetries to only 
lepton mass suppressed terms. Incidentally, we also present two procedures to introduce lepton 
mass corrections inside the observables in the distribution, an exact one preserving the structure 
of the massless case as much as possible and an approximated one that minimizes the error of 
neglecting lepton masses. The results presented here aim at reducing substantially two main 
sources of systematic errors of the experimental data namely the S-wave pollution and lepton 
mass corrections. The final goal is to provide strategies to enhance the sensitivity to any possible 
signal of New Physics affecting the clean observables Pi,2,3, P'a 5 6- 

In Section I we recall the definitions of the known observables, introduce new ones and we 
also present a so called "massless-improved limit" . In Section II we show how to extract the Pi 
observables in an exact way using folded distributions including all lepton mass corrections with 
zero S-wave pollution. Also a less clean approximate procedure is discussed using uniangular 
distributions. In appendices A and B we define the coefficients of the distribution coming 
from the companion scalar decay B — > Kq(—> i^7r)/ + /~ and provide more examples of folded 
distributions. 
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I. Definitions: Observables, full distribution and massless-improved limit 



The angular distribution that describes the four-body decay B — > K*(—> Ktt)1 + 1 including the 
S-wave pollution from the companion decay B — > Kq(-^- Kir)l + l~ is [6] l9l ITB] 



d 4 T 



dq 2 cfcos 9k cfcos #/ 1 



9 
32^ 



J±s sin 2 Ok + J\ c cos 2 9k + (J2s sin 2 Ok + J2c cos 2 9k) cos 29 i 



+ J3 sin 2 9k sin 2 9\ cos 2<f> + J4 sin 2$^ sin 29\ cos </> + J5 sin 2(9 ^ sin 0/ cos 

+( J6s sin 2 6* a- + Jq c cos 2 cos 0; + J7 sin 2#fc- sin 9\ sin + Js sin 29k sin 20j sin < 



+ Jg sin 2 sin 2 9\ sin 20 



(1) 



where 



S = — 



1 

4tt 



Jf a + Jf fe cos 9 k + (^|a + J 2b cos ^if ) cos + J a sin sin 29 £ cos ( 
+ J5 sin 9k sin cos + J-j sin 0# sin Be sin </> + Jg sin Bk sin 2#^ sin < 



and 



X = j dm K jBW K *(m 2 K7T 



2 \|2 



(2) 



(3) 



being a correction introduced in [16J to take into account the width of the resonance (see [16] 
for precise definition and details on this function). This correction factorizes from the standard 
Jj defined in [51 [71 [9] and as, we will show, it will always cancel in the Pi observables. The Jj 
containing the interference terms coming from the B — > Kq(-^- Ktt)1 + 1~ decay computed in [16] 
are defined in Appendix A. 

The full differential decay distribution is then 



dT 



full 



dq 2 



dT K * dT K> 
dq 2 dq 2 



with 



dTK* 1 <ITk* 
^ 2 = 4 (3«7ic + 6Ji s — J2 C — 2J 2s ) X, ^ 2 



(4) 



(5) 



Here we will not consider scalar contributions to facilitate the comparison with [16] where 
they were also neglected. Notice that the distribution including lepton masses (but neglecting 
scalars Jg = 0) contains 11 Jj coefficients (only 10 are independent PHI 15]) plus 8 Jj coefficients 
coming from the interference with the S-wave polluting process. A generalization of the results 
presented here to include scalars is straightforward. We take the same inputs as in [28] and 
all observables should be understood as J — > J + J (obviously the same applies for expressions 
written in terms of transversity amplitudes). 
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The first important point is the exact definition of the transverse asymmetries in terms of the 
coefficients of the distribution. In [16| the transverse asymmetries are defined with an unusual 
normalization 



7(2) _ 4J 3(9 ) A^)( n 2 \ - 4J 9(g 



3Ju(q 2 ) ~ hM 2 ) ' 3J ls (g 2 ) - J 2s (g 2 ) ' 

4"V)= o, r^rW ^ 

3Ji s {q 2 ) - J 2s [q 2 ) 

This choice of normalization in |16j is driven by the aim of extracting it from the Ox uniangular 
distribution (see below). These definitions albeit correct have two main disadvantages: 

• The expression of the transverse asymmetries given in Eq.([6]) once expressed in terms of 



transversity amplitudes involves lepton mass terms via the coefficien 



3xpr 
f 



Ju 



+ 



2\ 



\ A L\2 + l A Ll2 + lA R l 2 + lA R l 2 



4m 



2 



q 2 



+ ^ Re ( A L A R* + A L A R*\ (7) 



where = y 1 — 4m 2 /q 2 . This implies that their definition is not invariant under changes 
of the lepton mass (£ = e,n). 

• Normalization requires to measure two coefficients of the distribution. 

For these two reasons we prefer to stick to the definition of the Pi observables as given in 
|15j . These definitions in terms of transversity amplitudes do not dependent on any lepton mass 
term: 

\Ai\ 2 - \Af(\ 2 + (L o R) j 



Pi 



\ A L\2 + \ A L\2 + ( L ^ R) 2 J 2s ' 



Re(^Mjf -AfAf* ) Jgs 

P 2 — i ,r, 9 , . 7TT5 , TZ ^7 — Pi 



\A L ± \ 2 + \A\ | 2 + {L^R) ™8J 2s ' 

Im(^Mf-^f) ^ 
3 |^i| 2 + |Ajf | 2 + (L^R) 4J 2s ( } 

These observables were completed in [28] with 

y/2Be(A$*Af + A R A R *) j 4 



P' 

1 A 



A^\ 2 + \A^\ 2 + (L o (|A£| 2 + (L o i?)) 



\J—J 2c J 2 s 



p , = V2Re(A%*Aj-A R Af*) = (3, J 5 

! ' A L ± \ 2 + \A^\ 2 + (L^R))(\A^\ 2 + (L^R)) 2 



V2Im(A%*Af + A R Af*) ft J? 

^6 - ; o / j 7 I 9 ) 



^| 2 + |Ajf| 2 + {L^R) \ (\A%\ 2 + (L o i?)) 



2 V—'hc-hs 



1 From now on we will distinguish between two types of lepton mass dependences: lepton mass terms that will 
refer to terms like the last term in Eq.0 and mass dependence via the prefactor fit- 
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where the f3g prefactor included in P2,5',6' is there to ensure that the expression in terms of 
transversity amplitudes does not contain any lepton mass dependence. P\ is the transverse 
asymmetry A^p proposed in [6] and P2 and P3 correspond to the A^f 1 and A^ m \ respectively, 
normalized as in [14] and not like in Eq. (|6|) . Notice also that the normalization of the transverse 
asymmetries in Eqs.([8|) requires to measure one single coefficient of the distribution, namely, 
J2 S , free from any lepton mass term. 

It is very convenient to express the coefficients of the angular distribution directly in terms 
of observables. In Ref. [TS] it was presented the most complete expression of the coefficients of 
the distribution including lepton masses and scalars in terms of the full basis of observables. 
Later on, in [28] a very compact parametrization was found for the massless case without scalar 
contributions. Here we present a generalization of Ref. [28J to include lepton masses with a 
similar structure to the massless case at the price of trading the original basis [15] 

The exact parametrization of the coefficients in the massive case is ther@ 

3 * 1 dT K , - 1 «2 F 1 dTK " 



4 1 X dq 2 ' Zs 4 He 1 X dq 2 ' 

__ p 1 dr K * 2 ~i dv K * 

1 dTx* T n2rt ~ 1 dTx* 



J7 



PiHy/Vi ^ dJ & . ( 10 ) 



x 

and the redundant but interesting coefficient Jg is 



1_ 2 _, f—~ 1 dT K * 



Js = - 2 ^F T F L --^ (11) 



2 Notice that in the first coefficient Ji B the dependence does not appear explicitly (but it is included inside 
Ft). The reason is that parametrized in this way and only in this coefficient the mass dependence piece of fit will 
help later on to reduce the contribution from the lepton mass term. 
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where the new observables Ftl, Ft,l are given by 
"(2 + #) 



Ft 

F l 
F T 

F L 
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\A t \'' + 2Re(A^Af 
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3N 
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16J 2s 

/3|iV 



4 
iV 



^||^| 2 + K| 2 ] 



4J 2c 



(12) 



and the normalization is 



N 



(3 + /3|) ( \A L ± \ 2 + |A[| 2 + |4f + |A«| 2 + |^| 2 + \A R \ 2 ) + 



12m 



2 r 



A t \' z + 2Re(A^* + A^A^ + AjVA 



3Jic + 6Ji s - J; 



2c 



2J- 



2,s 



4 dT K * 
X da 2 



At this point several remarks are in order: 

I. Pi and (/ 2 -dependent observables once expressed in terms of transversity amplitudes do 
not dependent on rri£, on the contrary Ft^l and Ft,l do depend. 

II. Eqs. (|10p and Eq. fjlip are exact concerning lepton masses (no scalars included). In all the 
folded distributions discussed in next section we will always use this exact parametrization. 

III. They are constructed to resemble Eqs.(14) of [28; for the massless case. The main difference 
with the massless case is that here there are two longitudinal polarization fractions and 
two transverse polarization fractions and the sums of transverse and longitudinal fractions 
differ slightly from one: 

F L + F T = 1 + 6 F T + F L = l + 5 

In the exact massless limit F L \ me= Q = F L \ me=0 , F T \m e =o = F T \ mi =o and 5, 8 ->■ 0. A 
numerical expression for 5(s) is provided in Table 1. 5(s) will be re-expressed in terms of 
other quantities below. 

IV. The fit dependence in the coefficients is kept explicit (except for J\ s ). This choice will be 
fundamental when approaching the massless limit to reduce the systematic associated to 
this limit. 

We will explore now in a bit more of detail the massless limit. It is clear that in the exact 
massive case one should strictly fit the data to four observables Fl,t, Fl,t apart from the six 
clean Pj. The ratios Rt = Ft I Ft and Rl = Ft /Fl are indeed a direct measurement of the 
clean observables M\ and M 2 defined in |15| : 

R T = \ (2 + (1 + m x )pf) R L = l + M 2 
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Figure 1: (left) ratio R T = F T /F T always below 1.2% in SM, (middle) ratio R L = F L /F L with 
an error range between 1% to less than 4% in the SM (right) 5 represents the deviation from 
one of Fl+Ft- 



The deviation from one of these ratios in the SM in the range 1 to 6 GeV 2 are shown in 
Fig.l. They are below 1.2% for Rt and below 4% for Rl- Moreover, given the tight relation 
of these ratios with M\i and the relatively small deviations found in these observables in the 
New Physics scenarios analyzed in [15], a similar range of deviations is expected also for New 
Physics. 

Finally, since it is very costly in statistics to fit for the four observables Ft,l, Ft,l one could 
take the massless limit (mi — > everywhere) to reduce the number of observables to fit. However, 
the exact massless limit could induce errors of order 4% (around q 2 = 1 GeV 2 ) for instance if 
Fl = Fl is used (see Rl in Fig.l). For this reason we prefer to define an intermediate limit, 
that we call "massless-improved limit". This consists in keeping the dependence of Eqs, (|10p 
and reexpressing Ft and Fl in terms of Ft in the following way 



F T = F T + (Rt - 1)F T = F T + L\ 

F L = 1 - F T + 5+(R L -l)(l- F T ) + 0((R L - 1)5) = 1 - F T + L 2 



(13) 
(14) 



where 0((Rl — 1)5) is much below the permille level and 5(s) = L1 + L2. A numerical expression 
for Li ; 2 in the SM can be found in Table 1. In this improved limit the contribution to Flt 
from lepton mass terms (Li,!^) m t ne SM is always below 1.3% in all g 2 -region from 1 to 6 
GeV 2 . Notice that the Pi are not directly affected by this limit since they are always defined 
by the same Ft,l, but they are only indirectly affected via the propagation of the error in the 
determination of Ft,l- The reason for the sizeable reduction from 4% to near 1.3% in the case 
of Fl is due to the partial cancellation between 5 and (Rl — 1)(1 — Ft) in the SM (see Fig.l). 
This means that the error done by expressing the coefficients of the exact angular distribution 
Eqs. ffTOj) and Eq. ffTTT) only in terms of Ft,l (keeping the fii dependence but discarding the 
corrections) is below 1.3%. 

Concerning the integrated observables there are two important observations when including 
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lepton masses. First notice that 



dq 2 \ dq 

and 



~ dT K * ( ~ dT K * 



FtFl 'Esi = ( Jp~F l ^ 



dq 2 \ v dq 2 

which means that these products are the same for the massive and massless case. Second 
important remark is that three of the g 2 -dependent Pj observables include a fig in their expression 
in terms of the Ji coefficients and this can be problematic when defining integrated observables. 
The reason being that what is measured are the Jj coefficients and not products like figJ%. 
Consequently, a small redefinition of three of these observables is required when including lepton 
masses in the integrated observables. This means that the integrated observables in the massive 
case naturally split in two categories: 

• First set: Pi, P3 and P4 their g 2 -expression in terms of spin amplitudes and Jj coefficients 
do not contain neither lepton mass terms nor fig prefactors. 

ra e ^0 . _ kin J 3 _ fbin^e P ^T^T 



bin 



dr 



2 Ibin J ^ f bin Pe F T^ 



, pmt^O _ fun J 9 _ kin $\ ' P 3 F T 

^ ^3 >bin ~ . r j - _ a2i ,_ dT 



4 Ibin J ^ kn^ F T- 



K* 



dq 2 



< p ,ra^ ^ = kin J 4 = I Mn $ F i V ^F L ^ 

\f~ kin J 2c kin J 2« \j kin Pe F T^T kin Pl F L^T 

For these observables P™^° = P u P™^° = P 3 and = P' 4 . Naturally, only when 

expressing them directly in terms of observables the fig prefactor arises. 

Second set: P2, P5, Pg also their g 2 -expression in terms of spin amplitudes do not contain 
any lepton mass terms or fig prefactors. However their q 2 — differential expression in terms 
of Ji do contain a fig prefactor, which means that its corresponding massive definition for 
the integrated observables requires a redefinition with a prefactor 1/ fig: 



^ r 2 ^b 



2 ^ knPl F T%- 

k m h _ JunfliWMVwiP 



p/m^O _ Jbin J 5 _ J bin ^1 ^ J dq 2 

>bin — 



Ibin J 2c kin J 2s \/ kin $1 F ? Ibin $1 F L 

ir. 



p/mffu _ J bin _ - - aq- 

2 V ~ Am J 2c /bin J 2s ^ / w „ fig F T kin $1 F L 

where P™^° = P 2 /&, ^ m ^° = P^/A and P^° = 
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5 


dim 


GeV° 


GeV" 2 


GeV" 4 


GeV" 6 


GeV" 8 


GeV" 10 


c 

Oi 


-0.0498563 


+0.0502874 


n aooooa'7 

-0.0228207 


+0.00540332 


-0.000644906 


+0.000030558 


Rn 


+1.02522 


-0.0074101 


-0.0142605 


+0.00717305 


-0.00120443 


+0.0000692328 


Rhi 


+1.09905 


-0.09191 


+0.041595 


-0.00990168 


+0.00118789 


-0.0000565414 


L\ i 


+0.0201856 


-0.0273282 


+0.0135174 


-0.00331905 


+0.000403496 


-0.0000193199 


L2i 


-0.00598948 


+0.0255013 


-0.0149247 


+0.00389817 


-0.000488089 


+0.0000237576 



Table 1: Coefficients of the polynomial approximation to the deviation parameters 5, Rt, Rl, Li 
and L 2 , such that 5 = J2i=o,5 S i s ^ r t = J2i=o,5 R Ti s ^ R L = Ei=o,5 R Lis\ L x = J2i=o,5 L i ^ 
and L 2 = Y,i=o,5 L 2iS l - 



The last important point here is to compare the SM predictions for the massive integrated 
observables defined here < p™ 1 ^ > with the massless ones computed in |28| . We will focus 
the comparison on the 1-bin region between 1 to 6 GeV 2 . For the first set due to the balance 
of 0j factors between numerator and denominator of Eqs. (|15|) together with the slow variation 
of this factor with s, one should expect a tiny difference between the massive and the massless 
prediction. Indeed we found that the massive prediction in the SM is only +0.8% larger than the 
massless one for Pi and and below permille level for P3. In the second set the comparison of 
the SM integrated observable in the range 1 to 6 GeV 2 in the massless and massive case is —0.6%, 
+1.4% and +0.8% for P2, P§ and Pg respectively. This implies that the massless predictions in 
the SM obtained in |28| are valid to an excellent approximation also in the massive case. 

In conclusion in this section we have proposed to use either the exact parametrization or 
an approximate parametrization so called "massless- improved" limit for the experimental fit 
(instead of the pure massless limit). This improved parametrization has the advantage of in- 
cluding the dominant contribution from lepton masses with an error when compared to the 
exact result that amounts only to a maximum of a 1.3% in the SM for the terms that contained 
Fl t i n the region between 1 to 6 GeV 2 while it is exact for the q 2 -dependent Pi observables 
(up to the propagation of the error in the determination of Ft£)- The practical consequence 
of this improved limit is to effectively reduce the basis of observables to one observable less 
{ ^ , Ft , Fl ; -Pl,2,3 ; -^1,5,6 } ( srnce we are taking M\ — > 0). If one neglects in addition 5 (keep- 
ing the /3i dependence) by using Fl ~ 1 — Ft an additional error associated to this second 
approximation (see Fig.l) should be added. In this case the previous basis gets reduced to 
eight observables. This step by step procedure will help to know precisely the size of the error 
associated to each approximation. Once this is done, next step is to compare with the theory 
predictions. We have shown here that the comparison of the exact massive predictions in the 
SM with the massless predictions [28] in the 1 to 6 GeV 2 bin for the integrated observables 
shows a deviation that in the worst case amounts to 1.4%. 
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II. Can one extract the observables Pl,2,3 and P4 5 6 free from S-wave pollution? 

In order to answer this question two aspects have to be analyzed. On the one hand, the intrinsic 
character of the observable with respect to the S-wave pollution and, on the other, the extraction 
procedure of the observable. 

Concerning the first aspect, one should distinguish the observables that does not suffer from 
S-wave pollution by construction from those that have an inherent structure that includes the 
pollution. It is clear that the products 



F T —^-, F L —^-, PiF T —^- P'\ F L F T 



dq 2 ' dq 2 ' dq 2 3 v dq 2 

but also Pj, Pj with i = 1, 2, 3 and j = 4, 5, 6 are all intrinsically free from this contamination. 
Notice that these products enter in the definition of integrated observables. Consequently the 
corresponding integrated observables are also intrinsically free from S-wave pollution. 

On the contrary, for instance, the Si observables [9] are normalized by the full differential 
decay distribution 

s < = '^r < 17 > 

dq 1 

and they will be directly affected by the S-wave pollution given that the full differential decay 
distribution is largely affected by this disease (see Eq.(|5])). Any other observable normalized by 
, like for example a longitudinal polarization fraction normalized by the full distribution 

(call it Fl)-, will be strongly S-wave polluted, however the product Fl -g^r is again free from 
pollution. 

Concerning the extraction procedure, in the following we will show that it is possible to 
extract these observables completely free from any S-wave pollution in the exact lepton mass case 
and, afterwards, we will present for completeness also an approximate solution using uniangular 
distributions. 



Exact solution including lepton mass corrections using folded distributions 

The crucial point here is to determine the exact procedure to extract the Pi observables. There 
are basically two choices: uniangular distributions or the full angular distribution. The major 
drawback of the former is that by integrating angles the distinct angular dependence of the 
unwanted J, L is lost and mixing (and pollution) with the interesting J, L is unavoidable. If, instead, 
the full angular distribution is used the major problem is that there is not yet enough statistics 
to have access to the full distribution. However, the clever idea presented in Ref.[lJ of using 
"folded" distributions, have the double advantage of increasing the statistics and focusing on 
a restricted set of angular coefficients. This allows to extract information from this mode in a 
more selective way than traditional uniangular distributions and, in particular, it can be used, 
as we will show, to isolate the polluting Jj terms. In this note we will exploit this trick to show 
that using folded distributions one can measure all Pi and Pj free from any S-wave pollution 
thanks to their distinct angular dependence. 

In Ref. [1], the identification of 4> «-)■ 4> + tt (when <f> < 0) has been used to produce a "folded" 
angle </> G [0, ir] in terms of which a (folded) differential decay rate dT = dT((j)) + dT{(j) — ir) is 
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obtained This folding exploits the angular symmetries of the distribution and reduce to a subset 
the coefficients entering each folded distribution. If we use the same folding procedure including 
the extra terms coming from the pollution channel B — > Kq(—> Ktt)1 + 1~ one finds the well 
known structure for the angular distribution together with an extra piece with a distinguishable 
angular dependence 



d 4 f 



dq 2 cfcos 8k cfcos 9\ < 



9 
16tt 



Jic cos 2 9k + Jis(l — cos 2 Ok) + J2c cos 2 9k{2 cos 2 9g — 1) 



where 



+<^2s(l — cos K )(2 cos 2 6t-l) + J 3 (l - cos 2 9 K )(1 — cos 2 9i) cos 20 



+ J6 S (1 — cos 2 9k) cos ^ + Jg(l — cos 2 0jc)(1 — cos 2 6£) sin I 



1 

2^ 



^la + ^16 COS 6 K + J£ a + J 2b cos °K (2 COS t - 1 



X + Wi (18) 



(19) 



The extra factor W\ concentrates all the S-wave pollution of this folded distribution. Its distinct 
angular dependence can be parametrized by 



W\ = aw + bw (2 cos 2 



1) + cw cos 9k + dw cos #k-(2 cos 2 



1) 



(20) 



Fitting for these four parameters {aw, bw-, cw, dw) one can measure and disentangle the S-wave 
pollutionQ But, more importantly, the different angular dependence associated to the relevant 
coefficients involved in the definition of the -Pi, 2, 3 permits an extraction of these observables 
in a way completely free from any S-wave pollution. It is clear from the definition of -Pi, 2, 3 
in Eq.© that the folded distribution provides access to J3, Jq s and Jg, numerators of Pi, P2 
and P3 respectively. But also their common denominator J2 S is identified with the coefficient of 
(1 — cos 2 #if)(2cos 2 9i — 1). The normalization of A^, and Aj?^ presented in |16] requires 
to measure in addition J\ c , unnecessary for the Pi,2,3- 

This is more easily seen writing down Eq. (|18p in terms of observables using Eqs. (|10p . The 
following expression generalizes Eq.(l) of the LHCb note [1] by including lepton mass corrections 
and the S-wave pollution in a compact form: 



d A r 



dq 2 cfcos 9k cfcos 9i 
1 



9 
16tt 



F L cos 2 K + jF T (l ~ cos 2 9 K ) - $F L cos 2 9 K (2 cos 2 t - 1) 



+~$F T (1 - cos 2 6» x )(2cos 2 1 



1) + 2^l P l^(l - cos 2 K ){1 - cos 2 0t) cos 2<j) 



+2f3iP 2 F T (l - cos 2 K ) cos 9 e - /3|P 3 P T (1 - cos 2 9 K )(l - cos 2 9 e ) sin 24> 



^+w> ( 21 ) 



3 Other examples of singled folded distributions exhibiting dependencies on -Pi, 2, 3 but with the angle Ok instead 
can be found in Appendix B. 

4 Notice that if lepton mass suppressed terms from the polluting decay B — > KqII are disregarded aw = —bw 
and cw ~ —dw 
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or in the massless-improved limit: 
d 4 f 9 



dq 2 dcos 9 K cfcos 9 t d(f> 16vr 
1 



(1 - F T ) cos 2 9 k + jF T (l - cos 2 9 K ) - p$F L cos 2 9 K (2 cos 2 6> £ - 1) 



+^/3|F T (l - cos 2 ^)(2cos 2 ^ - 1) + -/3|PiF T (l - cos 2 0*-)(l - cos 2 9 t ) cos 20 



+2(3 £ P 2 F T (1 - cos 2 6>x) cos € - $P 3 F T {1 - cos 2 0*r)(l - cos 2 0<) sin 20 



+ (22, 



where the dependence is kept, Fl is traded by Ft and the L\ )2 subleading terms have been 
dropped off. The error of this improved limit is below 1.3% as discussed in the previous section. 

There is also an alternative that can be automatically applied to all the folded distributions 
discussed here and it is to use dT d f g f instead of , i-e, 



dq 2 



d 4 T 
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dq 2 dcos 9 K dcos 9 X # 167T 
1 



F L cos 2 K + ^Ft(1 ~ cos 2 9 K ) - 0\F L cos 2 9 k (2 cos 2 i 



1 



+-^f F T (1 - cos 2 K ) (2 cos 2 t -l) + -0\P x F T (\ - cos 2 9 K )(1 - cos 2 9 e ) cos 20 



+2/3 e P 2 F T (l - cos 2 Ok) cos 9 t - /3 2 P 3 F T (1 - cos 2 9 K ){1 - cos 2 9 t ) sin 20 
where all quantities with a bar are defined by Y = CY with 

C 



< ^ + W l (23) 



dq 2 



dTic* 



d£full 
dq 2 

Notice that, as explained in the previous section, no S-wave pollution is introduced in this way in 
the relevant q 2 or integrated P{ observables. Here the improved-massless limit becomes defined 
by 



Ft = CFt = C Ft + (Rt ~ l)C-fr = CFt CL\ 



(24) 



where the relative error induced by dropping CL\ in CFt is obviously the same as before, 
below 1.2%. However for Fl two possibilities are now open depending on the size of C. First 
possibility, if \C — 1| is comparable in size to CL 2 taking into account that C — 1 is negative 
while CL 2 is positive a partial cancellation can take place and the best approach is to use 



Fl = CF L = 1 - CF T + CL 2 + (C - 1) 



(25) 



In practice this require |C — 1| to be below 4% to be a good choice (given that L 2 is below 1%). 
The second possibility is to take 



F L = CF L = CF L + (R L - 1)CF L 



(26) 
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In this case the relative error — 1 is independent of C exactly as for Fr, however it rises 
up to 4% at q 2 = 1 GeV 2 as can be seen in Fig. 1. This choice is safer because the error 
does not depend on C but the high remaining error (4%) implies that the advantage of the 
massless-improved limit in front of the massless case is lost. 

At this point it is evident the power of this folding technique, in the massive case the full 
distribution would depend on 11 Jj coefficients (only 10 of them are independent) plus 8 Jj 
coefficients from the companion decay. Instead, the folded distribution depends only on seven 
Jj coefficients (or seven observables Ft, Fl, Ft, Fl, P\, P%, P3 ) and four Jj (a\y , b\y , c w , dw) ■ 
If the improved-massless limit is used the number of observables reduces to six: dTx* /dq 2 , Ft, 
Fl, Pi, Pi, P3 plus the four from J. 

Notice also that the prefactor X = f dm 2 Kn \BWK*(m K7T )\ 2 cancels out exactly in the Pj 
given their definition in terms of ratios of coefficients. 

One can also exploit further this technique and find several doubled-folded distribution show- 
ing different dependencies on the -Pi, 2, 3 observables and exhibiting common/different S-wave 
pollutions. For instance, two possible foldings of the angles <p an d ®t are 

IV. Identifying <p 4> + 71 when <p < and 0g o 9g — \ when 9g > | with (f> G [0, tt] and 
9 1 £ [0, 7r/2] such that the folded differential distribution is 

df = dT(4>, h, 0k) + dT{4>, h + 1,9 K ) + dF(4> - tt, § t , 9 K ) + dT(0 - tt, 9 t + \,9 K ) 



corresponding to 
d A f 



dq 2 dcos 9k dcos 9i dip 32tt 



4F L cos 2 9 K + 3Pt sin 2 9 K + F T sin 2 9 K {PjPi cos 2p> 

Tk 
dq 2 



+4f3 E P 2 (cos 9 1 - sin § e ) - 2$P 3 sin 20)1 + W A (27) 



where 



1 

W 4 = - 

TT 



J\a + Jib cos Ok 



V. Identifying <\> o <p + tt when <\> < and 9# O tt — 9g when 9g > | with ip G [0, tt] and 
9 1 G [0,7r/2] and the folded distribution 

df = dT{4>, 9 e , K ) + dY{4>, ^ - 9e, &k) + dT(4> - tt, 9 e , 9 K ) + dr(0 - tt, tt - 9 e , 9 K ) 
corresponding to 

d 4 f ! 



dq 2 dcos 9k dcos 9\ d<p 327r 



Acos' 9 K (F L - (3fF L cos 29 1 ) + sin 2 9 K (3F T + f3fF T cos 29 e )+ 

dq 2 



+2/3|sin z ^P T (Picos20- 2P 3 sin 20) sin^ £ ) — ^ + W 5 (28) 



where W 5 = 2W 1 (<p,9 e ,9 K )- 
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Pi 


P2 


^3 






n 


(1)-(13) 


(l)-(4), (6)-(10),(12),(13) 


(1)-(7),(H)-(12) 


(10) 


(10)-(13) 


(11).(12) 



Table 2: Sensitivities to each P observable from the set of 13 folding distributions discussed. 



Other double foldings with the angles <fi and 9k that offer a different dependence on the observ- 
ables Pi, 2,3 can be found in Appendix B. 

Finally, we show how the pollution can be also disentangled from the P' A 5 6 observables using 
this folding technique. Here we present different possibilities: 

X. Identifying the angle <p <H> — <p when <p < with <p G [0, 7r] and the differential folded 
distribution is now dT = dT((f>, 6g, 9k) + dT(—(f>, 9g, 9k) where 

— * = [4 cos 2 9 K (F L - $F L cos 29 1 ) + (3F T + 0jP T cos 26 e ) sin 2 9 K + 

dq 2 dcos9 K dcos6id(f) 647r L 

+2F T {0jP l cos 20 sin 2 6>k sin 2 9 t + 4/^P 2 cos 0* sin 2 9 K )+ 



+2\] ' F L F T (l3jP' A sin 29 K sin 20* + 2/3^ sin 29 K 



sm (A> cos ( 



^ + .^(29) 



with 



W w = — 



2vr 



Jfa + ^2a cos 29 e + cos 9k (Ju + <^2b cos 20^) + cos (f)( J5 + 2 J4 cos 0^) sin 0^ sin 0^ 



XI. Identifying the angle 0^ -H- 7r — 0^ when 0£ > | with 0^ G [0, 7r/2] and the folded distribution 
(if = dT(cp, §i, 9 K ) + rfT(^, 7T - 0*, 0*:) is then 



d 4 T 



4 cos 2 X (F L - 0jF L cos 20^) + (3P T + $F T cos 20*) sin 2 9 K 



dq 2 dcos 9k dcos 9 X d<f> 647T 

+20jF T (P l cos 2<p sin 2 X sin 2 0^ - 2P 3 sin 2(f) sin 2 e sin 2 0^) + 

+4/3 £ \J F L F T (P^ cos - P@ sin p>) sin 20^ sin ( 

where 



ii 



2vr 



Jf a + J 2 a cos 20^ + cos 9 k {J ib + «^2fe cos 2^) + (J5 cos <p + Jj sin 0) sin 0^ sin 0^ 



XII. Identifying 9g -H- 0^ — \ when 0^ > | with 0^ G [0,7r/2] and the folded distribution 
(if = dT((f>, § e , 9 K ) + dT((f>, h + *,0k) is 

d 4 f 9 



dq 2 dcos 9k dcos B\ dip 32ir 



(4F L + 3P T + (4P L - 3P T ) cos 29 K )+ 



+F T (-0fPi cos 2(f) sin 2 9 K + 2(3 e P 2 sin 2 0^(cos - sin § t ) - /3^P 3 sin 20 sin 2 9 K ) 



+/W F L F T (P^ cos (/)- Pq sin (p) sin 20 x (sin 9 e + cos € ) 
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where 



W 



12 



47T 



2Jf a + 2Ji b cos 6*x + ( J5 cos 4> + J7 sin <?!>) sin 6k (cos 0£ + sin #^ 



One last example of two-folded angle distribution sensitive to P§ comes from identifying <p -R- 
when (ft < and 6g 6e — ^ when 0g > S. The corresponding combination is 



7T 



<ir = (ir(^>, e K ) + dr(<t,,e, + -,e K ) + dr(-<p, e,, e K ) + dr(- 



where 



d 4 r 



tig 2 cfcos 6*x dcos 61 < 



with 



32vr 



-(AF L + 3F r + (4F L - 3F T ) cos 29 K )+ 



W 13 = — 
2ir 



+2f3 e \J F L F T P' h cos <j> sin 26^(sin Q t + cos 6g)+ 
+F T sin 2 6 K {fifPi cos 20 + 4/3 £ P 2 (cos ^ - sin 6 t )) 



2Jf a + 2 Jf b cos 9k + </s cos </> sin 9k (cos &£ + sin 6g) 



dV 



K* 



dq 2 



+ W 13 



Of course, using the same procedure also the integrated observables can be extracted cleanly 
from the folded distributions in a way completely free from any S-wave pollution. 



Approximate solution using uniangular distributions 

In Ref[16j the 'unusual' normalization factor 3 Ji s — J2s of the transverse asymmetries is obtained 
from the uniangular distribution 



d 2 T 



dq 2 d cos 6k 



a e K (q 2 ) + be K (q 2 ) cos 6 K + cq k (q 2 ) cos 2 6 K 



(30) 



where the coefficient functions are easily obtained integrating the angular distribution Eq.([T]) 
over the angles (j) an d 9g . In particular, in our notation one finds 



ae K (q 2 ) = J C la-\j2a + ^ @JlB -J2s)X = Jl a - X -J^ a +^{ 9f V - i 2 / j 



3_ 

32 



dTK* 
~dq T 



3 

32 



9F T - /3 2 F T ) (1 + A(q 2 )) 



dT K * 
dq 2 



(31) 



which corresponds exactly to Eq.(26) of Ref.|16| and defines in our notation the function A(q 2 ) 
introduced in Ref.[T6j. This function A(g 2 ) is the responsible for the deviation up to 23% 
around q 2 = 2 GeV 2 . Indeed it is not surprising this huge effect. It is immediate to see from 
the definition of A(q 2 ) in [16] in terms of transversity amplitudes that it is proportional to 
\M'q\ 2 + \M'^\ 2 (see Ref.[H>] for definitions) which is the source of this pollution and large 
deviation. Notice that with the procedure described in [16] this pollution is not suppressed by 
any m 2 /q 2 prefactor. 
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However, if one uses instead the observables Pi, P 2 and P3 denned as in Eq.([8|) one can design 
a procedure to extract their normalization J 2s using uniangular distributions with a lepton mass 
suppressed S-wave pollution. In the region between 1 to 6 GeV 2 this amounts to a much small 
pollution than the one found using A(q 2 ) coming from unsuppressed lepton mass terms. 

We should emphasize that using uniangular distributions is not the best choice, since as we 
have shown folded distributions allow for an exact extraction of the Pi free from any pollution. 
Still if one insist in using uniangular distributions a procedure to obtain this normalization is 
easily constructed using the distributions on 0£ and (p. The idea is to combine the two uniangular 
distributions 



d 2 T 



dq 2 d cos 0£ 
d 2 T 



dq 2 



ae e + be e cos Ot + c 9e cos 2 9 e 



+ cos <p + sin <fi + d^ cos 20 + sin 2(f) 



(32) 



(33) 



where the coefficients can be trivially obtained by integrating the corresponding angles. Two 
possible combinations relevant for our purposes are: 



a e + 37r 



ae, + 



3 + /3f 
l+l 2 



C0. = K 2 



(34) 



where 



#1 



K, 



1 



8(3+$: 

3 l K1 



3 + 



SPi (3F T + F T + 2(F L -F L ) 



dT K * 
dq 2 



+ 16(/3 2 J 1 c a + J 2 c J 



(35) 



In turn they can be written as 

3 



Ki 



XJ 2s + O 



K 2 = 3XJ 2s + O 



ra£ 
q 2 



which shows that both combinations allow the extraction of the normalization coefficient J 2s 
with an error given by terms of order 0{m 2 /q 2 ). 

The influence of the 0{m 2 /q 2 ) term is moderate in the region of interest between 1 and 6 
GeV 2 (below 1 GeV 2 more difficult problems and uncertainties arises). 



Conclusions 

We have shown that the Pi observables are not only intrinsically clean observables from S-wave 
pollution point of view, but also that they can be extracted using folded distributions in a 
way completely free from any S-wave contamination coming from the companion decay B — > 
KqI + 1~ , including all lepton mass corrections. The same conclusion applies to the corresponding 
integrated observables. 
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We have also defined a massless-improved limit for the differential distribution, an interme- 
diate stage between the massive and the massless case, that effectively reduces the number of 
observables at a relatively low cost: an error below 1.3% in the SM in the region between 1 to 6 
GeV 2 . Also a comparison of the SM prediction for the integrated observables < P{ >g 2 =i— 6GeV 2 
in the massive and massless case shows a discrepancy below 1.4% in the SM, implying that the 
SM massless predictions for the integrated observables given in[28] are an excellent approxima- 
tion. 

Both results are important experimentally to get a better control or even remove the sys- 
tematics associated to S-wave pollution and lepton masses. 
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Appendix A. J i: definitions 

The definition of the Jj derive directly from the full angular distribution given in |16| 



J c la =l?(q 2 ) / |B%.K,)rdmL 



JL = W) / \BW K *(m K7r )\'dm 2 K7r 



J{ h = 2^/3 J Re 
J| 6 = 2V3 J Re 
J 4 = 2^3 [ Re 
J 5 = 2V3 I Re 



J 7 = 2V3 j Im 
J 8 = 2V3 I Im 



l?BW K * (rn 2 Kn )BW^(m 2 K7: 
l c 2 "(q 2 )BW K * (m 2 Kw )BW^ (m 2 Kn ) 
Xl{q 2 )BW Kl { m 2 K ^BW^rn 2 Kl ,) 
Z'l{q 2 )BW K *{m 2 K ^)BW^{m 2 Kn ) 

Z'l{q 2 )BW K *{m 2 Kv )BWlS™ 2 K,) 
l^q 2 )BW K *(m 2 K7I )BW^(m 2 Kn ) 



dm 2 Kn 



dm 2 



Kit 



dm 2 KlT 



dm 



dm 2 



K 71 



dm 2 



Kt> 



(36) 



We refer the reader to [16] for definitions of Z; and Breight-Wigner resonances BWK*(m KlT ), 
BW K *(m 2 K7T ). See also [29]. 
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Appendix B. More folded distributions 
Observables Pi,2,3 

• Single folded distributions: The angle Ok can be folded in two different ways. 

II. The first comes from the identification Ok ^ it — Ok when Ok > § with Ok £ [0, 7r/2] 
implying the combination dT = dT(cj), Oi, Ok) + dT((j), 0£, ir — Ok) corresponding to 



d 4 T 



9 



dq 2 dcos K dcos 0\ dcj) 647r 



4 cos 2 K (F L - 0fF L cos 20 e ) + sin 2 K {3F T + 0fF T cos 26 e )+ 



+2F T (j3jP 1 cos 2(f> sin 2 0^ sin 2 e + 4j3 e P 2 cos e sin 2 0r- 



P3 sin 2<p sin 2 0^ sin 2 0£ 



dTK* 
~dq T 



+ w 2 



(37) 



where 



Wo = — 



2vr 



Jf a + J 2a cos 20^ + (cos 0( J5 + 2 J4 cos 0^) + ( J7 + 2 Jg cos 0^) sin 0) sin Ok sin ^ 



III. The second is Ok 0x — f when 0^- > ^ with 0/<- G [0,7r/2] and the combination is 
df = dT(<j), e , K ) + dT(</>, 0,, 0^ + |) given by 



d 4 r 



9 



<ig 2 dcos a- dcos 0^ d0 1287T 



4F L + 3F r + (-4F L + Ft)0£ cos 2^+ 



+2F T ($P 1 cos 20 sin 2 ^ + 4/3^ P 2 cos £ - 2/3fP 3 sin 20 sin 2 £ 



dr 



A'* 



dg 2 



+ W 3 (38) 



where 



2(J 1 c a + Jf a cos 20^) + ( Jf b + J| fe cos 20<?) (cos 0# - sin K ) + (cos 0( J 5 + 2 J 4 cos 1 



+(J7 + 2Js cos 0^) sin 0)(cos Ok + sin 0a-) sin( 



Double folded distributions: 



(39) 



VI. Identifying o + 7r when < and 6>x o Ok — § when 0x > § with G [0, 7r] and 
0A" G [0, 7r/2] with a combined distribution 



df = dT(0, e , K ) + dT(4>, e , 0* + |) + dT(4> - vr, 9 e , K ) + dr(0 - vr, e , K + |) 



given by 



d 4 T 



9 



dg 2 dcos K dcos X dcj) 647T 

+20fF T {P l cos 20 - 2P 3 sin 20) sin 2 



4F L + 3 Ft + (-4F L + F T )pj cos 20 e + 8F T (3 E P 2 cos e + 

dTK* 



dq 2 



+ W 6 



(40) 
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where 



Wa = 



2tt 



2(«/i C a + JL cos 26 f ) + ( J{ h + J c 2h cos 20 1 ) (cos 0^ - sin 9 K ) 



VII. Identifying o <p + 7r when < and 9k ^ tt — 9k when 0^; > | with G [0, it] and 
0_a; G [0,7r/2] and the combined distribution 



dT = dT(</>, e , 9 K ) + dT{4>, e , ir-9 K ) + dT(cf> - vr, 9 t , 9 K ) + dT(<j> - it, 9 e , it - 9 K ) 



where 



4 cos 2 9 K (F L - f3fF L cos 26 i) + sin 2 9 K (3F T + ^F T cos 29 £ )+ 



dq 2 dcos 9k cfcos 9i d(f) 32ir . 

+2 sin 2 6k Ft (pj{P\ cos 2cf> - 2P 3 sin 2<p) sin 2 6 e + 4/3^P 2 cos 9tj 



^ + W 7 (41) 



with 



VF 7 = 



l r~ 



7T 



^l C a + ^aCOs2^ 



VIII. Identifying <p <H> —0 when <p < and 9k ^ tt — 9k when 0k > f with </> G [0, 7r] and 
0_a; G [0,7r/2] with the combined distribution 



df = dT(cf>, e e , Ok) + dT{</>, 9 e , it-9 k ) + dT(-<f>, e , K ) + dT(-$, 9 e , vr - 9 K ) 



where 



d 4 t 



4 cos 2 K (F L - $F L cos 29 e ) + {3F T + $F T cos 26 1 ) sin 2 6 K 

dTK* 



dq 2 dcos 9k dcos 0\ dip 327r 

+2F T {l3jP 1 cos 2<f> sin 2 K sin 2 9 t + 4/^P 2 cos ^ sin 2 K ) 

with 



dg 2 



+ V^8 (42) 



1 



7T 



Jf a + J| a cos 20^ + cos 0( J5 + 2 J4 cos 9e) sin 0^ sin ( 



IX. Identifying <p ^ —cp when < and 9k 9k — § when 0^ > § with <p G [0, 7r] and 
9k S [0, 7r/2] with the combined distribution 



= dT($, e , K ) + dT(j>, e , 9 K + J) + dT(~4>, 9 e , K ) + dT(-$, 9 t , 9 K + J) 



where 



d 4 T 



4F L + 3F T + $(-4F L + F T ) cos 20 £ 



tig 2 cfcos 9 K dcos 9\ d<p 647r 

+2F T ($P 1 cos 2^ sin 2 £ + 4/3 e P 2 cos 1 



(43) 



with 



Wg = — 

2tt 



2J{ a + 2 J| a cos 20£ + ( Jf 6 + J| 6 cos 20 £ )(cos 0^- - sin X ) 



+ cos 0( J5 + 2 J4 cos 9p) (cos 0x + sin 0^) sin 0^ 
Finally, notice that the list is not exhaustive and other foldings are also possible. 



(44) 
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